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Abstract 

Origin of hydrodynamic turbulence in rotating shear flows is investigated. The particular em- 
phasis is the flows whose angular velocity decreases but specific angular momentum increases with 
increasing radial coordinate. Such flows are Rayleigh stable, but must be turbulent in order to 
explain observed data. Such a mismatch between the linear theory and observations/experiments is 
more severe when any hydromagnetic/magnetohydrodynamic instability and then the correspond- 
ing turbulence therein is ruled out. The present work explores the effect of stochastic noise on 
such hydrodynamic flows. We essentially concentrate on a small section of such a flow which is 
nothing but a plane shear flow supplemented by the Coriolis effect. This also mimics a small section 
of an astrophysical accretion disk. It is found that such stochastically driven flows exhibit large 
temporal and spatial correlations of perturbation velocities, and hence large energy dissipations of 
perturbation, which presumably generate instability. A range of angular velocity (f2) profiles of 
background flow, starting from that of constant specific angular momentum (A = Or 2 ; r being the 
radial coordinate) to that of constant circular velocity (iu = Or), is explored. However, all the 
background angular velocities exhibit identical growth and roughness exponents of perturbations, 
revealing a unique universality class for the stochastically forced hydrodynamics of rotating shear 
flows. This work, to the best of our knowledge, is the first attempt to understand origin of insta- 
bility and turbulence in the three-dimensional Rayleigh stable rotating shear flows by introducing 
additive noise to the underlying linearized governing equations. This has important implications 
to resolve the turbulence problem in astrophysical hydrodynamic flows such as accretion disks. 
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INTRODUCTION 



There are certain rotating shear flows which are Rayleigh stable but experimen- 
tal/observational data argue them to be turbulent. In absence of magnetic coupling, such 
flows are stable under linear perturbation. What drives their instability and then turbulence 
in absence of linearly unstable modes? In order to understand a plausible route to hydro- 
dynamic turbulence in such flows, a series of papers was published by different independent 
groups [ill. Based on "bypass mechanism" [6-8], those papers showed that such flows ex- 
hibit large transient energy growth of linear perturbation. It was argued earlier that indeed 
the large transient energy growth is responsible for subcritical transition to turbulence in 
plane Couette flow and plane Poiseuille flow, as seen in laboratory [9j. In such flows the 
transient growth increases with increasing Reynolds number (R e ). For example, in plane 
Couette flow, the maximum transient growth scales as R\. However, the situation changes 
drastically when rotational effects come into the picture. The Coriolis force is the main 



epid 



culprit behind this change. It was shown earlier [4J, |5|, [10] that rate of the increase of tran- 
sient growth in two-dimensional (insignificant vertical scale height) linearized rotating shear 
flows with Keplerian angular velocity (fl) profile (Q ~ r~ 3//2 ; r being the radial coordinate of 
the flow) decreases much compared to that in plane Couette flows, when maximum growth 
scales as R 2 J^ . Furthermore, in three dimensions (with finite height), such transient growth 
in Keplerian flows is insignificant, independent of viscosity, to generate any instability and 
turbulence. But the three-dimensional Keplerian flow is an important natural flow, which 



exists in several astrophysical contexts. Interestingly, while some authors [11], based on 
the prototype experiment consisting of fluid confined between two corotating cylinders, re- 
ported the Rayleigh stable rotating shear flows, similar to that of a Keplerian disk, laminar, 



121 ] found it turbulent. However, the results from direct numerical simulations 



some other 

\\\ argue that laboratory flows are indeed hydrodynamically unstable and should become 



turbulent at low Reynolds numbers. But the last author 



131 ] and the present authors 



also showed/argued that the experiments are compromised by undesired effects due to the 
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finite length of the cylinders. Note that some (effective) nonlinear theory 15|, ll6| revealing 
large energy growth of perturbation (even of finite amplitude) argues for plausible routes to 
turbulence in three dimensions. Hence, the real challenge is to uncover the mystery of mis- 
match between theory and observation under linear theory, in absence of a finite amplitude 
of perturbation. 

As argued earlier 10] that in three dimensions, one requires to invoke extra physics, e.g. 
external or self-generated noise etc., to reveal large energy growth or even instability in 
the system. It is needless to mention that the real flows must be three-dimensional, e.g. 
an astrophysical accretion disk must have a thickness, however small may be, and hence 
rotation plays a significant role to determine the dynamics of the flow. The aim of the 
present work is to investigate the amplification of linear perturbations which may eventually 
lead to hydrodynamic instability and then plausible turbulence in certain rotating shear 
flows supplemented by stochastic fluctuations in three dimensions. 

There is a long term association of growing, unstable modes generated by perturbed flows 
with statistical physics. The flow fields perturbed either at the boundary or through exter- 
nal forcing have been shown to produce emergent instabilities [3], altering the scaling 



structure of the systems. It was shown by Forster, Nelson and Stephen [19] that in the long 
time, large distance asymptotic limit, large non-equilibrium fluctuations in unbounded flows 
decouple to stabilize the flow. The authors also studied the critical dimension beyond which 
fluctuations get redundant. However below that dimension, all the Brownian oscillations 
oerturb the flow in the equilibrium limit. In a follow-up work by De Dominicis and Martin 
20J, a further generalization of such 'forced Navier-Stokes flows' was provided incorporat- 
ing a range of possibilities for the perturbing random forces which is correlated with the 



Kolmogorov spectrum 
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22] to redefine the scaling laws in the presence of such stochastic 



components. Later on, the focus was on studying topical variations in the transport prop- 
erties of fluids under stochastic perturbations [23] . These works were the precursor to the 
Kardar-Parisi-Zhang (KPZ) model 24J which studied a differential variant of the Burger 



field to i 
trum 
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lumine fundamental hydrodynamic instabilities related to the intermittency spec- 
, |26]. The KPZ model also sheds light on varieties of non-equilibrium phenomena, 
including fluid flows, bacterial growth, paper burning front, etc. (see, {27]). All such results 
were later encapsulated in another illustrative work by Medina, Kardar, Hwa and Zhang 
where the effects of varying spatio-temporal correlations were studied on the Burger 
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equation (which is a one- dimensional variant of the Navier-Stokes equation). It also showed 
for the first time that the statistical field theory could be a legitimate consort of fluid me- 
chanics such that local instabilities seen in fluid flows could often be identified as non-trivial 
infra-red divergences. All the studies, however, were restricted mostly to trivial/periodic 
and open/absorbing boundary conditions. 



Chattopadhyay & Bhattacharjee [29] applied the randomly stirred model of Navier-Stokes 
flows constrained by geometry in a flow involving a three-dimensional bounded incompress- 
ible fluid perturbed by boundary layer shear. In that work, they considered two relatively 
shearing flat plates with a layer of fluid sandwiched in between. The effect of the con- 
strained flow is implemented through the stochasticity structure function, which ensures 
that the symmetry is violated through noise autocorrelation function into the problem. 
Note that the geometry could also be changed over to rotating co-axial cylinders {soj for 
externally impressed shearing forces, essentially bearing the same qualitative imprint. This 
problem, known as boundary layer turbulence, has recently been revisited experimentally in 
laboratory, by generating shear in the flow through a raft of blowing air (mI- A numerical 
simulation governing large Eddy also reveals same experimentally observed non-inertial ef- 
fects [32j. Note in a related context that astrophysical observed data indeed arg ued for the 



signature of chaos in rotating shear flows, more precisely accretion flow s |33l . l34j|. which fur 



35| and neutron stars |36| due to turbulent 



ther supports infall of matter towards black holes 
viscosity. 

In the present study, we implement the ideas of statistical physics discussed above to rotat- 
ing shear flows in order to obtain the correlation energy growth of fluctuation/perturbation 
and underlying scaling properties. In the next section, we discuss the equations describing 
the stochastically forced perturbed flows which are to be solved for the present purpose. 
Subsequently, in §3 we investigate the temporal and spatial correlations of perturbation in 
detail, in order to understand the plausible instability in the flows. Finally, we summarize 
the results with conclusions in §4. 

II. EQUATIONS DESCRIBING FORCED PERTURBED ACCRETION FLOW 

As our specific interest lies in hydrodynamic instability and turbulence in presence of 
stochastic noise, we straight away recall the Orr-Sommerfeld and Squire equations in pres- 
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ence of the Coriolis force 



10], but supplemented by stochastic noise. They are estab- 



lished by eliminating pressure from the linearized Navier-Stokes equation with background 
plane shear (0, —x, 0) neglecting magnetic coupling in presence of angular velocity Q ~ r~ q in 
a small section of the incompressible flow [44|]. Hence, the underlying equations are nothing 
but the linearized set of hydrodynamic equations combining with the equation of continuity 
in local Cartesian coordinates (x,y,z) at an arbitrary time t (see, e.g., 0, 10 1 for detailed 
description of the choice of coordinate in a small section) given by 

d d\„ 2 2d( 1 _ 4 



d d \ _ du 2 du 1 „, _ , 

where u, ( represent respectively x— components of velocity and vorticity perturbation, R e is 
the Reynolds number, 7] 12 are the components of noise arising in the linearized system due to 
stochastic perturbation such that < r)i(x, t)rjj(x', t') >= Di(x)5 3 (x—x')5(t — t')5ij. Note that 
equations and (j2J) describe perturbations in a small section of an accretion flow which 
could be expressed in local Cartesian coordinates. The long time, large distance behavior 
of the velocity correlations are encapsulated in Di(x) which is a structure pioneered by 
Forster, Nelson and Stephen In the Fourier transformed space, the specific structure of 
the correlation function D^k) depends on the regime under consideration. It can be shown 



for all (non-linear) non-inertial flows 19|, |29( that D^{k) ~ l/k d , where d is the spatial 
dimension, without vertex correction and Di(k) ~ l/k d ~ a with a > 0, in presence of vertex 
correction. Note, however, that Di(x) is constant for white noise. 

We now focus onto the narrow gap limit which otherwise leads to an approximation x — > 
L, when L being the (narrow) size of the shearing box in the x-direction. The appropriate 



choice of a sharing box has already been described by the previous authors; see, e.g. |5l. |10|. 
which we do not repeat here again. In brief, in a local analysis one considers a small 
radially confined region of the accretion flow, while the azimuthal confinement is imposed 
by a periodic boundary condition. This effectively gives rise to a shearing Cartesian box 
for geometrically thin flows. We further choose a very small L so that we are essentially 
interested in an orbit. This is in the same spirit as of the celebrated work by Balbus & Hawley 



371 ] who initiated the existence of Magneto- Rotational- Instability (MRI) with axisymmetric 



perturbation. One of the motivations behind it is as follows. Presence of turbulence in 
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shearing astrophysical accretion disks is mandatory to generate viscosity therein. It is only 
the viscosity which helps in transporting the energy-momentum from an outer orbit to an 
inner orbit leading to the inspiral of matter, giving rise to a differential rotating (shearing) 
structure of accretion disks. As the aim is to address the origin of such viscosity (which can 
not be of molecular) for rotating shear flows, choosing x to be fixed justified. 

Hence, in the above discussed framework, we can resort to a Fourier series expansion of 
u, C and r)i as 



u{x,t) = J Zfoe^-^cPkdu), 
((x,t) = J C^e^-^^kdw, 
Vi (x,t) = [ TH^e'fa^tPkdu,, (3) 



where k and u respectively represent the wave- vector and angular frequency of perturbation, 
and substituting them into equations ([I]) and (jSJ) we obtain 

> )=M\ ^ | , (4) 

where 

M = 1 KA KA 1 ' (5) 

q 2 R e (k 2 — iR e u — iLR e k y 



Mi 

M 2 = 
M 3 = 



—k 2 q 2 (k 2 


— iR e u — iLR e k y ) 2 + 2(q — 


2)R\k 2 




2%qR\k z 




—k 2 q 2 (k 2 


— iR e u — iLR e k y ) 2 + 2(q — 


2)R 2 k 2 




i(q - 2)qR% 





-k 2 q 2 (k 2 - iR e UJ - iLReky) 2 + 2{q - 2)R 2 e k 2 ' 



j. k 2 q 2 R e (k 2 - iReto - iLR e k y ) 

4 ~ k 2 q 2 (k 2 - iR e u - iLReky) 2 - 2{q - 2)R 2 k 2 ' 1 ' 



when ffij^^', i = 1,2, are the components of noise in k — u space, k = ^Jk 2 + k 2 + k 2 . 

III. TWO-POINT CORRELATIONS OF PERTURBATION 

We now look at the spatio-temporal autocorrelations of the flow fields u and ( [27] . 
For the present purpose, the magnitudes and gradients (scalings) of these autocorrelation 



functions of perturbations would plausibly indicate noise induced instability which could 
lead to turbulence in rotating shear flows. 



A. Temporal correlation 

Assuming < Vikutf 2 ku >= ^> we obtain the temporal correlations of velocity, vorticity 
and gradient of velocity given below as 



< ufrt)u(x,t + r) >= C u (r) = J d kdu e— < u^u_ % ^ >, 

< ({x, t) C(x , t + r) >= C c (r) = J d 3 k dcu e*" < >, 

du(x,t) du(x,t + r) _ , . f , 3 , , _ ia;T ,2 
< — ^ — ^ " >= CduCO = / < ^u_ fc>a) >, 



(7) 

where r be the time at which the correlations to be measured. Now using equations (j4|), 
©, © and © 

C tt (r) = J d 3 kdcoe- luJT [2D iq A R 2 e (k A + (R e co + LR e k y f) 
+8D 2q 2 Rtkl]/[-A(-2 + q ) 2 Rtkt 
+AR 2 e q 2 {-2 + q)k 2 k 2 z (k 4 - (R e co + LR e k y ) 2 ) 

-kY(k 4 + (R e oo + LR e k y ff\. (8) 

We further consider the projected hyper-surface for which k x = k y = k z = fc/v3, without 
much loss of generality for the present purpose. As our one of the major interests is to 
determine the scaling laws, this restriction would not matter. However, this is equivalent to 
a special choice of initial perturbation which will affect the magnitude of the correlations. 
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We now perform the w-integration by computing the poles of the kernel given by 




2R 2 e (2-q) 2V2R e (2-q) _ LR e k 

3q 2 1 VSq f VS 

2R 2 e (2-g) _ kA _ 2V2R e (2-q) _ LR e k 
3g 2 1 V3q V3 h 

1 \ 2Rli2-q) ~ ( 2y/2R e {2-q) ^ LR e k, 



(9) 



By summing up the residues at the appropriate poles in presence of white noise, we 
evaluate the magnitude of frequency part of the integration given by 

[(bin — an) sin(a) cosh(o) + (am + bn) cos(a) sinh(fo)], 

where (10) 



a = y 2x + y/Ax 2 + y 2 , b = yJ x + ^/Ax^y 2 , 

X 2,(2-0),., 2 ,-. - k 2 



m = 2k\2-q), n = ^6(2=qjA. 



Hence 



fkm 

C u (t)=Aix dkk 2 RES Cu (q,R e ,r), (11) 
Jko 

with D X = D 2 = 1, 

where k = 2n/L ma _ x and k m = 27r/L min , when L = L max — L min being the size of the 
chosen small section of the flow in the radial direction, chosen to be 2 throughout for the 
present calculations. Similarly, C^{r) and Cd u { T ) can be easily obtained from equation (j7j). 
Note that the choice of constant Di and D 2 is due to our primary consideration of linear 
perturbation of the flow variables in presence of an additive noise. However, D\ and D 2 
need not be constant (viz. [19j) always. Later, we will discuss the cases with D ~ k~ d what 
effectively does not alter our results. 
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1. Asymptotic scaling laws 



At a large R e (— > oo) C u reduces to 




r, when q — >■ 2 or r -»■ 0, (12b) 



but ~ exp f — ^-r ) for large r. (12c) 

Hence, the correlation effectively appears independent of q at small and large r, at a large 
R e in particular. However, at an intermediate r, correlation diverges at q — > 2 — when the 
background specific angular momentum is conserved. Then with the decrease of q, due to 
appearance of fluctuation due to rotational/Coriolis effects with nonzero epicyclic frequency 
(see, 0,0] for a detailed discussion), correlation decreases (see Fig. [I] also). However, the 
correlation is effectively controlled by the dynamics of the flow, independent of viscosity. 
For completeness we also calculate the correlation for nonrotating flows, given by 

/■fem sinh ( ^ 

C u (T)nonrot = S^ 2 R e / dk (13a) 

Jk K 
T 

~ r when — << 1 (13b) 

Re 

~ R e when R e , r — > oo (13c) 
r 3 / 2 „/. [7\. . r 



it 



e 



erfi [kJ— |£™ when r ->■ oo or — » 1. (13d) 



Hence, unless r is very small at a finite R e , correlation depends on R e — quite an opposite 
trend of the rotating cases. Interestingly in absence of pressure, C u (r) nonrot ~ r 3 as r — > oo. 
Therefore, the rotation and pressure both create significant impact in the behavior in the 
flow such that the scaling behavior changes. 

Once w rms , %p and Cms (when A Tms 2 = | < [A{x, t + r) - A(x , t)] 2 > | ~ | < 
A(x, t) A(x, t + r) > |, where A = u, ||, C) are known, the evolution of perturbation energy 
growth £ = ![w rms 2 + ^( dU Qg s ) 2 + Cms 2 ] can be computed. However, by simple calculation of 
scaling law, it can be shown that each of the terms in S scales with r identically and hence 
the scaling behavior of £ is same as that of C u {r) apart from a constant factor. 
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FIG. 1: Temporal correlation as a function of time for R e = 10000, when q = 1 (solid line), 
1.5 (dashed line), 1.9 (dotted line), 1.9999 (dotdashed line) and flow is non-rotating (thick 

line). 
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FIG. 2: Temporal correlation as a function of time for q = 1.5, when R e = 10000 (solid 
line), 1000 (dashed line) and 100 (dotted line). 
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FIG. 3: Same as in Fig. [21 except for non-rotating flow. 
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2. Numerical solutions 



Figure [T] shows that with the decrease of q, temporal correlation of velocity perturbation 
decreases as the fluctuation due to the presence of Coriolis force increases. At small values 
of q (e.g. for 1, 1.5), however, the Coriolis fluctuation interacting with noise appears to be 
dominating in a way that the correlation reveals an oscillatory amplitude (this is more clear 
from Fig. [2]). Hence in a few rotation time the correlation exceeds that of large q, when it 
increases steadily in absence of dominant Coriolis effect at large q. However, in either of the 
regimes, temporal correlation of perturbation appears large enough in a few rotation time, 
which reveals nonlinearity in the system, presumably leading to instability. The correlation 
appears to be increasing steadily for the nonrotating case, in absence of Coriolis fluctuation, 
similar to that for q — y 2. Thus the interesting outcome here is the following. The Coriolis 
fluctuations along with a steady effect of noise lead to a larger C u {r) for a smaller q than 
that of a larger q in a few rotation time, as is also indicated by equation (fT2l . Note that 
an important physical quantity in this context is the rotation time of the flow. In a few 
rotation time, noise and the Coriolis fluctuation together bring in a large C u (t) and hence 
perturbation energy (prone to a smaller q, such as 1.5), which presumably could lead to a 
transition to turbulence. Therefore, the presence of noise makes the Coriolis effect active to 
generate instability, which were otherwise acting as the agents hindering the perturbation 
effects (see, e.g. [|3]), leading to a faster rate of perturbation energy growth into the 
system. 

Figure [2] shows that for q < 2, the correlation tends to become independent of R e . This 



is because with decreasing q, the epicyclic frequency {k = VL^j2{2 — q)) becomes larger 
bringing lots of fluctuation in the flow, which makes the flow limited by the dynamics itself 
but not by the viscosity. In absence of noise, perturbation energy growth (and hence C u (t)) 
is very small, as shown previously [4], for q < 2, particularly in three dimensions. Note, 
however, that at a small/moderate R e , epicyclic fluctuations work with noise in such a way 
that C u {t) at small r appears smaller than that of larger R e (see R e = 100 case in Fig. 
[2]) due to the harmonic nature of fluctuations. This effectively argues that the correlation 
saturates only with the increase of R e for q < 2. 

Figure [3] shows that C u (t) for a nonrotating flow at R e = 100 appears larger at large r 
compared to that at larger values of R e . This is due to the effect of noise, which would have 
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been suppressed at larger R e . This is also clearly understood from the analytical scaling laws 
given by equation ffT3l) . Note that, a very large R e along with a large r diminishes fluctuations 
in the flow leading to a linear perturbation growth with R e , as shown in equation (II 3d) . 

Comparing rotating and non-rotating cases, it can be further inferred that the noise 
suppresses the viscous effects in the rotating flow making it independent of R e in the regime 
of smaller r, as also shown by equation (112bl) . However, the optimum perturbation of a flow 
might reveal the correlation varying with R e , and the above inferences are valid for a fixed 
perturbation, considered for all the present cases. 

Note that in absence of noise, perturbation energy growth (and hence CAt)) in all the 
above rotating cases, essentially for q < 2, is very small, as shown previously [4(, particularly 
in three dimensions. However, the effects of noise bring a steady growth in the perturbation 
top of the Coriolis fluctuations. A remarkable feature in the scaling nature of all these corre- 
lations, atleast for asymptotic r (see, e.g., equations ( 112bj) and (I12cp ). is their independence 



of q (background angular velocity profile) — a trait identified in statistical physics literature 
as universality. All these flows, with identical overall perturbation energy growth exponents 



27l | in asymptotic limits of r (and also R e ), but generally with varying energy dissipation 



amplitudes, indicate a specific universality class for them. 



B. Spatial correlation 

Here also we assume < rfit w >= like temporal correlation and obtain spatial 
correlations of velocity, vorticity and gradient of velocity, given below as 

< u(x, t) u(x + f, t) >= S u (r) = J d 3 kduj e 1 ^^ < u% u u_%_ u >, 

< C(x, t)((x + f, t) >= S ( (r) = J d 3 k du e** < C fc >. >, 

du(x, t) du(x + f, t) 
dx dx 

where r be the position where the correlations to be measured. 

Now using equations (@J, fl3]), (jBJ) and (JT3|) . the velocity correlation perturbation S u (r) is 
explicitly given by 



< = >= SWr) = d^e'-V x < fi.fi. >, (14) 
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S u (r) = 2tt / dkk 2 d6 e ikrcose / du [2D iq A R 2 e 

Jk JO J 



'k 

4R 2 e q 2 (2 - q)k 2 k 2 z (k A - (R e u + LR e k y f) + A; 4 g 4 (fc 4 + 
(R e u + LR e k y ) 2 ) 2 ] 

rk m r 

= 2tt 2 / dk k 2 J (kr) / doo [2D iq 4 R 2 e (k 4 + (R e oo + LR e k y ) 2 ) 

Jk Q J 

+%D 2 q 2 R A e k 2 z ]/[A(2 - q) 2 R\k\ + 4R e V (2 - q)k 2 k 2 z {k A - (R e 

io + LR e k y ) 2 ) + k 4 q 4 (k 4 + (R e oo + LRk y ) 2 ) 2 ] , (15) 

where Jo(kr) is the zeroth-order Bessel function. Similarly, one can obtain S^(r) and Sd u {f) 
explicitly, when the poles in equation (fT5"|) are exactly the same as that given by equation 
([9]). Here also we stick to the simplifying assumption k x = k y = k z = fc/v3. 

Now following the same procedure as applied for the correlation of temporal velocity 
perturbation, we obtain the correlation of spatial velocity perturbation S u (r) given by 

rk m 

S u (r) = dkk 2 J (kr) RES 5u (g, R e ,r), (16) 

J k 

where KESs u (q, R e ,f) is the magnitude of the part arising from frequency integration. 



(k 4 + (R e u + LR e k y ) 2 ) + 8D 2 q 2 Rtk 2 z ]/[4(2 - q) 2 R\k\ + 



1. Asymptotic scaling laws 

At large r and R e (which is a natural/astrophysical scenario), the flows with all q reveal 
identical scaling of correlation with r given by 

kr). Hence, the 



S u (H e and r->oo) -> "^lyr ( 2 ~ <?) J dk exp(-kr) — ^ (17) 

which is independent of R e . Note that at a finite r, exp(— kr)/y/r — > Jq 
flow is limited by dynamics. This reveals identical roughness exponents 27J for all flows, 
indicating a single universality class, as was already predicted from the temporal correlations. 
Equation ( JTTj) also reveals a steadily damped instability. It further recovers the fact that the 
flow exhibits unbounded correlation as q — > 2, and hence unbounded growth of perturbation. 
At small r 

-km (2-q- 2k 2 ) log (f ) 



PKm 

S u (R e ^co, r^o) -> -Q (2 - <?r 3/2 / dk ■ 

Jkn 

(18) 



'ko K 
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FIG. 4: Spatial correlation as a function of radial coordinate for R e = 10000, when q = 1 
(solid line), 1.5 (dashed line), 1.9 (dotted line), 1.9999 (dotdashed line) and flow is 

non-rotating (thick line). 




FIG. 5: Spatial correlation as a function of radial coordinate for q = 1.5, when R e = 10000 
(solid line), 1000 (dashed line) and 100 (dotted line). 

which is again determined by rotational parameter q only, when \og(kr/2) becomes negative 
at a small r. Hence, except q = 2 (or q very very close to 2) and/or r = (or r very very 
close to 0), the correlation is finite. 

For completeness, we also describe scaling laws for nonrotating flows. At large r and R e 

q nonrot ^ ge^Vg exp(-fcr) 

°«(/? e andr^oo) ~^ fl / 2 aK ^9/2 ' ^ 

which scales with R e , unlike the rotating flows. Hence the flow is not limited by dynamics 
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FIG. 6: Same as in Fig. [5j except for non-rotating flow. 



at any r. However, at r — > the scaling behavior changes and is given by 



r k log- (Ml) 

Su^oo, r-0) f dk (20) 



which still scales linearly with i? e . 



5. Numerical solutions 



Figures HHS] show that the spatial correlations of perturbation decrease with the decrease 
of q from 2, when the epicyclic fluctuations arise into the flow. Moreover, the correlations 
for all R e tend to merge at large r, at a given q particularly, becoming independent of R e , as 
is also shown by the asymptotic result given by equation (TIT)) . However, at small r and finite 
R e , they depend on viscosity. This is due to the effect of noise, as is the case for temporal 
correlations. 

On the other hand, for non-rotating flows, the correlations are not limited by the dynam- 
ics of the flow, as shown by Fig. [61 rather are controlled by viscosity. This is also understood 
from equations (fl9|) and fT20l . Hence, the rotational effects constraint the correlation no- 
ticeably. 

It is generally seen that the correlations decrease with increasing r and decreasing q. Note 
their exponential decaying nature in equations f[T7|) and f TTTJ]) compared to the logarithmic 
variance in equations (ITS]) and f )20p . However, their value appears significant enough to reveal 
a steadily damped instability in the flow. The possibilities of steadily damped oscillatory 
instabilities in constrained Orr-Sommerfeld flows were already predicted in two dimensions 
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FIG. 7: Same as in Fig. SJ except for D 1 = D 2 ~ 



381 ] as well in the recirculating cavities between the jets produced behind a high solidity 



grid consisting of a perforated plate 



All the above correlations are computed for white noise. Figures [7] and [8] show the 
variations of spatial correlations of velocity perturbation for color noise respectively without 
and with vertex corrections in D^{k). Interestingly, while the magnitude of correlations 
decreases with respect to that for white noise (particularly for D\ = D 2 ~ l/^ 3 ), with the 
increase of q they appear large enough to govern instability. As an example, for noise with 
vertex correction, the Keplerian flow still reveals large enough correlation even at small 
radial coordinates. 



All the above large values of perturbation energy growth are indicative of pure hydrody- 
namic instability and plausible turbulent transport, which is one of the outcomes under this 
theory. 

C. Spectrum 

In order to compute the spectrum of flow, let us first integrate cu-part of equation (fl5l) 
and obtain 

j E u (k) dk = J d 3 kRES Su (q,R e ,k,D 1 ,D 2 ). (21) 
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FIG. 8: Same as in Fig. HI except for Di — D 2 ~ 1/k. 

Similarly, one can obtain J E^(k)dk and J Ed u {k)dk from and Sd u respectively. They 
would scale in a same way as of J E u (k)dk. Hence, the energy spectrum can generally be 
given by 

E u (k) = Ck 2 RES Su (Q,Re,k,D 1 ,D 2 ), (22) 

when C is a constant. Figures 191 and [TU1 show that the energy spectra for all the rotating flows 
with q < 2 have same spectral indices at large and small k. Moreover, the energy appears 
identical for all q at large k, which is due to the fact that a small length scale appears 
identical for all rotational profiles. At large scale, the spectral index is same as that of 
noise correlation, which implies that the flow is controlled by noise in that regime. However, 
interestingly, there is no change of slop in the nonrotating case, and the corresponding 
power-law index appears same as that of the large k rotating ones. This may be understood 
as the Coriolis force in the rotating flows distorts the system which results in affecting the 
large scale (and hence small k) behavior of the flow. Hence, the behavior of non-rotating 
flow appears same throughout without being affected by the rotational effects. Note also 
that the color noise increases the power-law indices throughout compared to that for white 
noise. 

IV. SUMMARY AND CONCLUSION 

We have attempted to address the origin of instability and then turbulence in rotating 
shear flows (more precisely a small section of it, which is a plane shear flow with the Coriolis 
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FIG. 9: Energy spectra for the velocity perturbation, when q = 1 (dashed line), 1.5 (dotted 

line) and 1.9 (dotdashed line) and flow is nonrotating (longdashed). The solid and 
thick-solid lines respectively indicate the nature of slop of the curves at small and large k 

for rotating flows. D 1 = D 2 = 1, R e = 10000. 




FIG. 10: Same as in Fig. El except for D X = D 2 ~ 1/k 3 . 

force). Our particular emphasis is the flows having angular velocity decreases but specific 
angular momentum increases with the radial coordinate, which are Rayleigh stable. The 
flows with such a kind of velocity profiles are often seen in astrophysics. As the molecular 
viscosity in astrophysical accretion disks is negligible, any transport of matter therein would 
arise through turbulence. However, the Rayleigh stable flows, particularly in absence of 
magnetic coupling, could not reveal any unstable mode which could serve as a seed of tur- 
bulence. Therefore, essentially we have addressed the plausible origin of pure hydrodynamic 
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turbulence in rotating shear flows of the kind mentioned above. This is particularly mean- 
ingful for charge natural flows like accretion disks around quiescent cataclysmic variables, 
in protoplanetary and star-forming disks, and the outer region of disks in active galactic 
nuclei etc. where flows are cold and of low ionization and effectively neutral in charge. Note 
that whether a flow is magnetically arrested or not, the magnetic instability is there or 
not, hydrodynamic effects always exist. Hence, the relative strengths of hydrodynamics and 
hydromagnetics in the time scale of interest determines the actual source of instability. 

We have shown, based on the theory of statistical physics, that stochastically forced 
linearized rotating shear flows in a narrow gap limit reveal a large correlation energy growth 
of perturbation in presence of noise. While the correlations are directly proportional to 
the Reynolds number for nonrotating flows, they are independent or weakly dependent of 
Reynolds number for rotating flows. Hence the later flows are mostly controlled by the 
dynamics and noise. While correlations of perturbation decrease as the flow deviates from 
the type of constant specific angular momentum to of the Keplerian and then to of constant 
circular velocity, still they appear large enough to trigger nonlinear effects and instability. 

While earlier the origin of hydrodynamic instability in rotating shear flows of said kind was 
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addressed based on "bypass method" in two dimensions 0, Q] , it was shown that 
physics, e.g. elliptical vortex effects, has to be invoked for such instabilities to occur in three 
dimensions. However, later theories showing instabilities are effectively based on nonlinear 
theory. Therefore, the present work addresses the large three-dimensional hydrodynamic 
energy growth of linear perturbation for the first time, in such rotating flows, to best of 
our knowledge, which presumably leads to instability and subsequent turbulence in such 
flows. Only requirement here is to have the presence of stochastic noise in the system, which 
is quite obvious in natural flows like astrophysical accretion disks. Indeed, to best of our 
knowledge, no one so far has addressed the effects of stochastic noise in such rotating flows. 

Interestingly, the flows with velocity profile Q ~ r~ q , with q < 2, exhibiting identical 
overall growth exponents (but may not be identical energy dissipation amplitudes), indicate 
a specific universality class. At large r, they also reveal identical roughness exponents. 
Thus the properties of temporal and spatial correlations together indicate the Rayleigh 
stable rotating shear flows to follow a single universality class. 

Let us end by clarifying suitability of astrophysical application of our results which is 
based on incompressible assumption. If the wavelength of the velocity perturbations is 
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much shorter than the sound horizon for the time of interest (e.g. infall time of matter), 
then the density perturbations (i.e. sound waves) reach equilibrium early on, exhibiting 
effectively uniform density during the timescale of interest. For a geometrically thin as- 
trophysical accretion disk around a gravitating mass, the vertical half-thickness of the disk 
(h) is comparable to the sound horizon corresponding to one disk rotation time 40] such 
that h ~ c s /Q. Therefore, as described in previous work , [16|, |4l|, |42j, for processes 



taking longer than one rotation time, wavelengths shorter than the disk thickness can be 
approximately treated as incompressible. A detailed astrophysical application of the present 
work has been reported elsewhere 4^ |. 
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Appendix: Magnetic version of the set of Orr-Sommerfeld and Squire equations 

Let us now establish the linearized Navier-Stokes equation in presence of background 
plane shear (0, —x, 0) without neglecting magnetic coupling with a background magnetic 
field (0,i?i,l), when Bi being a constant, in presence of angular velocity Q ~ r~ q in a 
small section of the incompressible flow. Hence, the underlying equations are nothing but 
the linearized set of hydromagnetic equations including the equations of induction in local 
Cartesian coordinates (see, e.g. 
in a small section) given by 

d d \ 2v dptat 1 ( '„ dB x dB x \ 1 , 



10| for detailed description of the choice of coordinate 



dt X dy) U q ^ dx An \ 1 dy dz ) R e ^ U} 
d ^d\_ (2 1 V , dp to t 1 ( r, 9B y dB y \ 1 2 



9 _ 9 \ . . , dptot 1 f n 9B Z dB z \ 1 ^ 2 



m- x 9y) w + ^-TA B ^ + Hi (25) 
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(26) 



dB v dv „ dv dB x dB z „ 1 ~ „ . . 

when the vectors for velocity and magnetic field perturbations are (it, v, w) and (B x , B y , B z ) 
respectively, R e and R m are the hydrodynamic and magnetic Reynolds numbers respectively, 
p to t is the total pressure perturbation (including that due to the magnetic field). Above 
equations are supplemented by the conditions for incompressibility and no magnetic charge, 
given by 

du dv dw . . 

dx dy dz 

dB x dB v dB z , , 

ox ay dz 

Now by taking partial derivatives with respect to x, y, z respectively to both the sides of 
equations ( |23|) . ( 12^1) . ( |25i) and adding them up we obtain 

2 2 f dv du\ f n du 



^ ^ tot q \dx dy) dy ^ ^ 

Now taking V 2 to equation (|23|) . using equation ( )3T|) and defining x-component of vorticity 
( = dw/dy — dv/dz, we obtain 

d d\ ^ 2 2d( 1 ^ n 3 , 1^ 4 



Also by taking partial derivatives with respect to x and y respectively to both the sides of 
equations fl2^|) . (123]) respectively, subtracting one from other and defining ( B = dB z /dy — 
dBy/dz, we obtain 

d <9V / 2\ du 1 / d <9V 1 2a 

Furthermore, by taking partial derivatives with respect to y and z respectively to both the 
sides of equations (1281) and (|27|) and subtracting one from other, we obtain 

8 .. 8 \ , 8( in 8(.8B, , 1„, 



a-^J & = ai + B 'aJ + lF + ^ v& - < 34 > 
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The equations (152]) . (!33|) . ( 126|) and (134")) describe the set of magnetized version of the Orr- 
Sommerfeld and Squire equations in presence of the Coriolis force, for the background mag- 
netic field described above and linear shear. To best of our knowledge, this set of equations 
has not been reported in the existing literature so far. 
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